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Abstract 

We define a modification of the Erdos-Renyi random graph process which 
can be regarded as the mean field frozen percolation process. We describe 
the behavior of the process using differential equations and investigate their 
solutions in order to show the self-organized critical and extremum properties 
of the critical frozen percolation model. We prove two limit theorems about 
the distribution of the size of the component of a typical frozen vertex. 



1 Statements 

The frozen percolation process on a binary tree was defined by D. J. Aldous in 111: 
it is a modification of the percolation process which makes the following informal 
description mathematically rigorous: we only occupy an edge if both end- vertices 
are in a finite cluster. The self-organized critical property of this model manifests 
in the fact that for f > 5, which is the critical time of the corresponding percolation 
process, a typical finite cluster has the distribution of a critical percolation cluster. 

I. Benjamini and O. Schramm showed that it is impossible to define a similar 
modification of the percolation process on 1?. An explanation of this non-existence 
result can be found in Section 3. of |2j|. 

First we give an informal description of the mean field frozen percolation pro- 
cess: It is a modification of the Erdos-Renyi random graph process: Initially we 
have a (not necessarily empty) graph on [A'^-mo(0)J vertices (one should think 
about A'^ as being lai^ge, but the initial mass mo(0) is fixed), and between every pos- 
sible pair of vertices, edges appear with rate Simultaneously lightnings strike 
vertices with rate X{t)iu{N) at time t and when a vertex is struck, the fire spreads 
along the edges and burns the connected component of that vertex: that subgraph 
is removed from the graph, including vertices. Thus the number of vertices of the 
random graph decreases with time. The expressions "burnt", "frozen", "deleted" 
and "removed" are treated as synonyms in the sequel. 

If l^^{t) denotes the number of vertices contained in components of size k 
in the random graph at time t, then the vector-valued stochastic process = 
{'U^ {t),l'2 {t),. ■ ■ ) also has the Markov property (the main advantage of the mean 
field model is that the graph structure of the connected components has no effect 
on the evolution of component sizes). We are interested in the model when 1 ^ A'^. 

Denoteby N = {1,2,...} andNo = {0,l,2,...}. 
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Definition 1.1. We fix nio(O) G M+. The mean field frozen percolation process on 
N vertices is a continuous time Markov process with state space 



= {V E N^ £ n < LA^- mo(0)J , Vfc ^ G No} 
k>i 

We define the coagulation and deletion operators 

+ ( {1^1, l^2,---,i^k-k,... ,1^1-1,. ..,l^k+i+k + l,.. .) ifk<l 

\ [1^1,1^2,. ..yk-2k,... ,'U2k + 2k,...) ifk = l ^' 

:=('J/i,...,V,-^,...) (2) 

Let X : M+ M+ be a positive continuous function and ^ : N ^ M+. The transition 
rates of the Markov process are 

+ N \ Ti-'^k-l^l ifk<l 

Mi^^^i") = MO-M^)-n- (4) 



Let v^{t) '■= -^^^ denote the mass of components of size k at time t. 

The mean field frozen percolation model is closely related to the mean field 
forest fire model (discussed in |[6l), the only difference in the definition of the 
Markov process is that in the case of the forest fire model, a burnt component 
of size k is replaced by k isolated vertices, so that the number of vertices in the 
random graph remains unchanged. The two models both have the self-organized 
critical property (and we believe that they are in the same universality class, which 
means that the theorems of this paper have analogous "forest fire" versions), but 
the corresponding partial differential equations have an explicit solution in the case 
of the frozen percolation model which enables us to say more about this model. 

CO 

V := {v = (vi)^^i : G M, > and £ v^t < 00} 

k=i 

\* ■={y: ■ve\,3K<+ooyk>K = 0} 

Definition 1.2. We consider a sequence of mean field frozen percolation processes 
with N ^ °°, but with the initial state 

„o).K,o).v»,o).^.„<,o),o,o.^.)-(^.^.....^.o,o,..^)ev 

and the lightning rate function X{t) fixed (independently ofN). Such a sequence is 
called 

• subcritical if^{N) = 1 
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• critical if jj <^ lj{N) < 1 

• alternating if fi{N) = jj. 

lfvk{^) = ^{k=\} ■'Wo(O) then the initial state is called monodisperse, otherwise it 
is polydisperse. 

We are going to describe the time evolution of the Umit object 

lim<(0=v,(0. (5) 

AT— >oo 

We introduce differential equations to characterize the limiting component size 
distributions vu{t) where G N and t G R+. They are modifications of the Smolu- 
chowski coagulation equation with multiplicative rate kernel: 

2 k—\ oo 

Ck{t) = :t £ ^ - • • Ck-i{t) - Ckit) £ / • c;(0) Flory's model 

^ 1=1 1=1 

(6) 

I k—\ oo 

Ck{t) = :^ £ ' ■ — 'Ck-i{t) — Ck{t) £ I ■ci{t) Stockmayer's model 
^ 1=1 1=1 



(7) 



If we let Vk{t) = k • Ck{t) then © becomes 



n{t) = :t £ v,(Ov;t-/(0 - ^-ViW- L v,(0) (8) 
^ l=\ k=\ 

We are going to use the formulation ([8]) rather than the classical 

The differential equations ([U) describe the time evolution of {vk{t))'^^i defined 
by ^ for the dynamical Erdos-Renyi random graph process (see |fTl|). If we only 
look at the evolution of the component size vector ^(?) in the dynamical Erdos- 
Renyi random graph model, we get the Marcus-Lushnikov process (see Q) with 
multiplicative kernel which is the ^^{N) = case of our model (no deletions, only 
coagulations). 

Definition 1.3. // {vk)1^^ =\£\ let 

mo := £ Vk mi := £ kvk mi := £ k'^Vk m^ := £ k^Vk 

k>l k>l k>l k>l 

Remark 1. Our definition of the moments m„ differs from the convention of the 
literature of the Smoluchowski equation by a shift of indices. 

If we define 

<{t) := £ vf (0 and 0^(0 := £ vf (0) - £ vf (0 = ^(^(0) -^5^(0 
i=\ i>\ i>\ 

(9) 
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then for all k the random function {t) is decreasing and {t) (the mass of burnt 
vertices) is increasing. 

It might happen (e.g. in the case of the Erdos-Renyi model) that 

6(0 := lim hm K(0 - (0) / lim lim K(0-wf(0) =0. 

In this case the mass missing from the small components is contained in a giant 
component of mass < Q{t). 

Definition 1.4. If\{t) is a solution of ([U), we define the gelation time by 

r« :=inf{? : mi(?) = +oo}. 

It is well-known from the theory of the Smoluchowski coagulation equation 
that an alternative characterisation of the gelation time is 

r« = inf{?:mo(?)<'Wo(0)}. 

For the solution of ([8]l the gelation time is = ;^7j|Qy, the mass of the giant 
component is 6(f) = mo(0) — mo{t). \{t) undergoes a phase transition: 

• For < f < the system is subcritical: 6(f) = and k^Vk{t) decay expo- 
nentially with k. 

• For < t the system is supercritical: 6(?) > and k i— > Vk{t) decay expo- 
nentially with k. Further on: f i-^ 6(f) is smooth and strictly increasing with 
limf^oo6(f) = mo(0). 

• Finally, utt = T^ the system is critical: 6(f) = and 

oo 

Y,Vk{T^)^K-^'^ as ^^oo. (10) 

k=K 

Our aim is to understand in similar terms the asymptotic behavior of the system 
when, beside the Erdos-Renyi coagulation mechanism, deletions due to lightnings 
also take place. 

Definition 1.5. We say that v(f ) = (v/t(f))^^j G V solves the general frozen perco- 
lation equation on [0, T] with initial condition v(0) G V*, a continuous nonnegative 
rate function X : M+ and control function <I> : M+ if 

yO<s<t<T < 0(0) < 0(5) < <I>(f) < mo(0) (11) 

and for all k = 1,2, .. . the equations 

vk{t)=vk{0)+ f'^'yvi{s)vk-,{s)-kvk{s){{mo{0)-^{s)) + X{s))ds (12) 
Jo 2 

and the inequality 

Vf 0< 6(f) :=mo(0)-mo(f) -<!>(?) (13) 

is satisfied. 
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It is easy to see by induction that the absolutely continuous functions vi (f) , V2 ) ) • • • 
are completely determined by ([121 ). the initial condition v(0) and the functions X 
and <I>. The only reason why we do not write 

Vk{t) = - £ vi{t)vk-i{t) - kvu{t) ((mo(0) - 0(0) + (14) 
^ i=\ 

instead of (fT2l) is that the increasing function <!>(?) might have jumps. 

There are three versions of the general frozen percolation equation correspond- 
ing to the three regimes on Definition 11.21 

• The subcritical system of integral equations are (fT2l) with the extra condi- 
tions Vf < Xinf <X{t) and 

<I>(0=mo(0)-mo(0- (15) 

That is 6(f) = by ([T3] ) (no giant components appear due to frequent light- 
nings) and the equations take on the form 

Vk{t) =Vk{0)+ -yvi{s)vk-i{s) -k-Vk{s)mo{s) -X{s)k-Vk{s)ds (16) 
Jo 2,tt 

The term —'k{s)k-Vk{s) indicates that in the subcritical regime even small 
components ai^e burnt with a rate proportional to their sizes and X{s). 

• The critical equations are (fT2l) with the extra conditions X(f) =0 and ([TSl l: 

Vk{t)=Vk{0)+ -yvi{s)vk-i{s)-k-Vk{s)mo{s)ds (17) 
Jo 2 f-[ 

X{t) = indicates that in the critical regime lightnings are not frequent 
enough to do any harm to small components, but (fTSl ) indicates that they 
are frequent enough to keep the mass of the giant component at zero. 

• Let = Tq <T^ < ■ ■ . he. a sequence with no accumulation points. Let 

M{t) := max{/ : t}' < t] (18) 
v(;) solves the alternating equations with burning times 7*, Tj,... if 

n{t) = \'Zvi{t)vk-i{t)-k.Vk{t)m^{T'^^^^) (19) 
^ l=\ 

Mind the difference between ([Hi and (fTTl ): in the case of the Erdos-Renyi model 
the small components are allowed to coagulate with the giant component (which is 
of size 6(?) = mo(0) — mo{t) by <!>(?) = and (fT3l)). but in the case of the frozen 



5 



percolation model the giant components are removed at the time of their birth. 
Using the terminology of the theory of Smoluchowski coagulation equations we 
might say that in the case of dH) the gel and the sol do react in the post-gel phase 
(Flory's model, Q), but in the case of ([IT] ) they do not react (Stockmayer's model, 
(IT])). Nevertheless, for t < the solutions of dS]) and ([TT] ) are identical since 
mo{t) = mo(0) in this regime. 

The intuitive meaning of ( fT9l ) is that giant components are removed from the 
system at the burning times. 

Thus OHl is ([HI) with 

e(0=mo(r^(,))-mo(0 (20) 

M{t) 

0(0 = mo(0) - mo(r^(,)) = mo(0) - mo{t) - Q{t) = £ e(r/) (21) 

7=1 

Both 9(f) and 0(f) are left-continuous functions of t. 

Note that in the case of the (sub)critical frozen percolation equations ( ([T6l l and 
(flT] )) the fact that <I>(f) is an increasing function automatically follows by (fTSl ): 

<t>{t) - <t>{s) = mo{s) - mo{t) = 

^ / --Y^vi{u)vk-i{u)+k-Vk{u)m(i{u)+X{u)-k-Vk{u)du = 

ft N oo 

lim / V V k-Vk{u)vi{u) +X{u) ■ k-Vk{u)du >0 



Theorem 1.1. 

• For any v(0) G V* and < < X{t) the equations (1161 ) have a unique 
solution. 

• For any v(0) € V* the equations (1171 ) /zave a unique solution. 

• For a«}' v(0) E V* anj sequence of burning times the equations ( 1191 ) 
have a unique solution. 

We prove this theorem in Section (3] 

Definition 1.6. The solution of the random alternating equations with rate func- 
tion X : M+ — > M+ is a \ -valued continuous-time Markov process: \{t) evolves 
deterministically, driven by the equations (1191 ). but the sequence of burning times 
7^,72, ... is random: 

Jim lp(f < r^(,)^i <t + dt\!Fr)= ^(06(0 (22) 
where ft the natural filtration generated by the process. 
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In plain words: a lightning strikes and burns the giant component with rate 
proportional to its size and 

Definition 1,7. 

'■= {{Wk)k=l '■ < Wl < W2 < • • • < +°°} 

■= { (wA:)"^j G W : 3K < +00 yk>K Wk = WK} 

Ifw G W denote by mo: = sup^w^:. 

We say that {{wk{-))'^^i ■•^i')) a frozen percolation evolution on [0,T\ with 
initial condition {wk{0))°^^i =w£'W*,or briefly 

if for all < t < T we have {wk{t))'7^i G W , for all k the functions Wk '■ [0, T] — > 
[0,mo(0)] are left-continuous and decreasing, : [0,T] — > [0,n2o(0)] is left con- 
tinuous and increasing with initial condition ^{0) = 0, moreover for all t <T we 
have (fT3] ). 

We define convergence on the space WyfjO, T]: 

(W(-))r=n^"(-))-(M-))r=n^(-)) 

as n^ 00 if for all k we have w'^{t) — > wtit) for all t which is a point of continuity 
ofwk and <!>"(?) (£>(t)for all t which is a point of continuity of^. 

With this topology the space 'JV„[0, T] is metrizable, complete and compact. 

From the frozen percolation process of Definition 11.11 one gets a random ele- 
ment of 'H'wiO, T] by Denote the probability measure on "W^wiO, T] correspond- 
ing to the process by P/y. 

It is easy to check that {{wk{-))'7^i ,^{-)) G '^yyjO,T] where Wk{t) =LLi v;(f) 
and \{t) is a solution of the general frozen percolation equation (fTTl) & (fT2l ) & (fT3l) . 

Theorem 1.2. We consider a sequence of frozen percolation processes (see Defi- 
nition [7T7]) with initial state \^{0) = v(0) G V* and X{t) positive and continuous. 
Define w^{t) and <^{t) as in (|9l). Denote the probability measure on 'H^yvjO, T] 
corresponding to the process by F^. 

Then P/v converges with respect to the weak convergence of probability mea- 
sures on the polish space Wy,[0, T] to a limiting measure P, which depends on the 
decay rate of fi{N) in the following way: 

• If ^{N) = 1 then P is concentrated on the unique solution of (1161 ) with rate 
function 

• ¥ jj i"!^) ^ 1 ^ concentrated on the unique solution of (I17I ). 

• If lLi{N) = then P is the law of the solution of the random alternating 
equation (see Definition \I.6i with rate function X{t). 
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We prove the iu{N) = 1 and the ^ <^ <^ 1 part of this theorem in Section|4l 
In fact, these proofs are almost identical to the corresponding convergence results 
of 161, but we present them here as well for the sake of completeness. 

We omit the proof of the iij{N) = part of Theorem ll.2l . but we believe that 
the methods introduced in SectionUl can be easily generalized for this case as well. 

If we formally substitute X{t) = into ([l6]l or T^^^^ = t into ^9^, we get (fT/l) . 
Rigorously: 

Theorem 1.3. Let (v"(?))"^i be a sequence of solutions of (1161) with the same 
initial condition v(0) G V* where X„{t) — > uniformly asn^oo. Then for all t and 
k lim„_>oo v^(?) =Vk{t) where \{t) is the solution of (1171) with the same initial data. 
lim„^oo<I>,j(f) = ^{t) uniformly on [0,oo). 

In plain words: if the rate of lightning is very small in the subcritical equations, 
then the solution is similar to that of the critical equation. We prove this theorem 
in Section m 

Theorem 1.4. Let (v"(?))"^j be a sequence of solutions of ( 1191) with the same 
initial condition v(0) where the sequence of burning times satisfy 



Then for all t and k lim„^co v^(?) =Vk{t) where \{t) is the solution of (1171 ) with the 
same initial data. lim„^oo = ^(t) uniformly on [0,oo). 

In plain words: if the burning times of the alternating equations are very fre- 
quent, then the solution is similar to that of the critical equation. We prove this 
theorem in Section |7] 

The solution of (fTTl ) has the self-organized critical property: for all < f it 
has the power-law decay of (ITOl) : 

Theorem 1.5. If\{t) is a solution of (1171) with initial condition v(0) G V*, then 
T^ = = Its ^crit{s)ds where (fcnt '■ [T^,+oo) —>■ is positive and con- 

tinuous, and for all t >T^ we have 



Definition 1.8. Letx*{t) := inf{;c : IX=\ Vk{t)(^ < +°°} 

The solutions of our equations have a remarkable rigidity property: 

Theorem 1.6. Ify{t) is the solution of ([T6l ) or ([T9l) and\{t) is the solution of ([TT] ) 
with the same initial condition, then for all t >T^ and k> I we have 



lim sup{7;^ (n) - T,\n)} = 0. 




(23) 



Vkit) = Vk{t)e 



-kx*{t) 
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The solution of ([TTl) with monodisperse initial condition is well-known (see 
e.g. HI) and explicit: 

Claim 1. If \{t) is the solution of (fTTl) with Vyt(O) = '^{k=i} •'Mo(O) then for t 



> 



= — Lt = — ^TvT we have 

mi(0) mo(0) 



1 k^-'^ , 

V.M = -—.-'■. (24) 

That is, for all < ? in the N ^ °° limit, the component size of a uniformly 
chosen (unburnt) vertex in the critical frozen percolation model has Borel distribu- 
tion, which is the same as that of a vertex in the Erdos-Renyi graph ait = T^. The 
Borel distribution ( (v^:(l))^^i in (|24l) ) is the distribution of the size of a critical 
Galton- Watson tree with POI{\) offspring distribution (see ifTl). 

The same self-similarity phenomenon can be observed in Aldous' frozen per- 
colation model (see HI) on the binary tree: for ? > 5, which is the critical time of 
the percolation process on the binary tree, a typical finite cluster has the distribution 
of a critical percolation cluster. 

The solutions started from a polydisperse initial state are asymptotically self- 
similar: 

Theorem 1.7. If \{t) is the solution of the critical equation (fTTl) with v(0) G V*, 
and vi (0) > then 

lim t-Vkit) = e and lim t-mo(t) = 1. (25) 

Theorems 11.51 . 11.61 and 11.71 are proved in Section fusing the method of 
Laplace transforms, which is classical for the Smoluchowski equation with multi- 
plicative kernel. The results (1251 ) and = Ck{t) k^^^^ (which is a variant of 
(I23] )) are already present in fSl, but we believe that our approach based on the no- 
tion of the critical core of \{t) (defined in Section |2]) gives new insight into these 
results about the solution of (ITTl) . 

In the frozen percolation model on the binary tree, components are frozen (i.e. 
removed from the system) when their size becomes infinite. The question may 
arise: 

What is the typical size of a frozen component in the mean field process of 
Definition \1.1V 

In order to precisely formulate this question recall (O and let 

^\[tut2lk):=^--\[te[tut2] ■ ^{t+)=!iU{t-)} ■ 
Thus <i>^ {[ti,t2\,k) is the mass of burnt components of size k from ti to ?2- We have 

£<I>^([ri,?2],^) =^>''(f2)-<J>''(?i) =:^i[h,t2]) 

k>\ 
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Thus p^[t\,t2] '■= %vf[ff ' k= 1,2,... is a random probability distribution for 
all N and ?i < ?2- 

Denote by | C^axi^) \ the size of the largest component at time t. 

Conjecture 1.1. If n{N) = N^^ in a critical sequence of frozen percolation pro- 
cesses (see Definitions \l.l\ and U.2i . where < a < 1, and if we define 



2a if a < i 
2 '/a > 3 



then for every <t we have 



log (E K(o)) 

hm ^\ = a (27) 

A?^»o log(Af) 

log(E(m^(0))-log(E(m7(0)) , 
lim ^ . . . ^ = P cx (28) 

log(A^) 

.in, ISiiEMlM) = (M, 

w^oc log(Af) 

Moreover for every v(0), < < f2 CC there exists a non-defective probability 
distribution function F : (0,°°) (0, 1), hm;t^o+^(-^) = 0, limx^ooF{x) = 1 such 
that for all x G M+ we have 

hm £l[^<x^P(«)].;.f[?i,f2] = F(x) (30) 

In plain words we might say that after gelation the typical component size of a 
frozen vertex and the size of the largest component is of order A/^P("). This conjec- 
ture is supported by heuristic arguments, computer simulations and Theorems 11.81 
and 1 1 .91 below. For < a < | the model is conjectured to behave similarly to the 
subcritical case described in Theorem 1 1.81 whereas for ^ < a < 1 it is conjectured 
to behave similarly to the alternating case described in Theorem 11.91 Note that 
P(^) = j and A'^S is the order of the size of the largest component in the critical 
Erdos-Renyi random graph. 

Theorem 1.8. If \^{t) is the solution of (1161) with rate function X{t) = X and 
\^{0) = v(0) G V* then there is a constant C that depends only on the initial data 
and T such that for all <X< I and ^^^^^-^ < t < T we have 

\Mt)-(?cntit)\<CX (31) 

where 

^^x{t)=^x{t)=lm\{t). (32) 
Moreover if we define the random variable Yx{t) to have distribution 
, , l-k-vl(t) k-vht) 
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then 

limP —Yy^{t)<x]= -^—e-^dy (33) 

In plain words: for any t > the distribution of the size-biased sample from 
the component-size distribution v^(f) rescaled by converges in distribution to 
a r(^, 1) distribution as X ^ 0. We prove this theorem in Section |7] 

The relevance of Theorem ! l.Sl to Conjecture [TTT]is the following: if we consider 
a sequence of subcritical frozen percolation models (see Definition ! 1.21 ) with = 
X then by Theorem ! 1.2! we get 



lim lim + = lim , ix — 

lim = ' = P (Yxit) = k) 

If we let X ^ then by dlB and we get m]{t)>iX^^ which is a "subcritical" 
version of (|27]), '4j\ = E {Yx{t)) x corresponds to p(a) = 2a in (EH), and ^ 
is a version of (l30l ). 

Theorem 1.9. Le? v^(0 denote the solution of the random alternating equations 
(see Definition U .6\ } with a constant rate function X{t) = X. 

Let 5{X) be a function satisfying X^2 <^ 5{X) <^ \ as X ^ °°. 

Recalling ([TSll and (l20l) let 

M(t) 

<i>,(?,x):= £e^(r,^)i[e^(r,^)>x] 

7=1 

be the random mass of frozen giants of size at least x. Then 



lim ^ s,,, ^ , , ^ ^ = / -^y^e-y dy (34) 

m probability. 

We prove this theorem in Section |7] 

The heuristic meaning of this theorem is the following: if we pick a vertex 
uniformly from all vertices that were frozen between t and t + h{X) and denote 
the mass of the giant component of that vertex by Zx{t), then the distribution of 



\\J -^^^"^xi^) converges to a size-biased Rayleigh distribution (see Definition !7.3! ) 

as X — > oo. Thus the typical mass of a frozen giant is of order X^ ^ , which suggests 
that if n{N) = jj- (that is a = 1 — £ in Conjecture !1.1! ) then the typical size of a 

frozen component is of order (N'^y^ -N = N^^^^, that is P(a) = ([34l) is the 
"alternating" version of (l30l ). 
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The critical frozen percolation model has an extremum property compared to 
the subcritical and alternating models (see Definition 11.21 ): if each burnt/frozen 
vertex produces profit at a rate ^$ per time unit after it has been frozen, but each 
lightning (even the ones hitting burnt vertices) costs jv.mo(o) ^' ^^^^ asymptotically 
(as N ^°°) the critical model is the best choice if we want to maximize our profit 
on [0, r]. We reformulate this extremum principle in terms of the differential equa- 
tions ([111), CH), CSl). 

The asymptotic value of our profit produced by burnt vertices as N ^ °° is 
Jq <t>{t)dt according to Theorem ! 1.21 The asymptotic cost of lightnings is Jq X{t)dt 
for the solution of (fT6] ). but it is zero for (fTTl ) and ([T9l ). since the price we have to 
pay for the lightnings vanishes in the case of critical and alternating models as 

N ^oo. 

Theorem 1.10. We fix v(0) G V*. Let y^" (t) denote the solution of (fTTl) with initial 
condition v(0) and let v*"^ denote the solution of (1161 ) with lightning rate function 
X{t) and the same initial condition. Then for any T > 

r ^'"\t)dt - r x{t) dt < r {t)dt - rodt (35) 

Jo JO JO Jo 

If\''^'{t) denotes the solution of (1191 ) with an arbitrary sequence of burning 
times and initial condition v(0) then 

r {t)dt < r {t)dt (36) 

Jo Jo 

Remark 2. Let T > = ^^^^ and £ > 0. For a suitable choice of'k{t) we have 

r ^"''\t)dt-{i-e) r x{t)dt> r ^"'•'{t)dt-{i-e) r m oi) 

Jo Jo Jo Jo 



For a suitable choice of burning times 

cT rT 



f ^""{t)dt + £^""{T)> f ^"-"{t)dt + e^"'''{T) (38) 
Jo Jo 

The idea that the critical forest fire model solves a variational problem is al- 
ready present in IS. 



2 Definitions, Transformations 

We consider a solution of the general frozen percolation equation (see Definition 

O). 

Denote the Laplace transform (generating function) of \{t) by 

oo 

Vit,x) := ^ Vkit)e-'" (39) 

k=l 
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for ;t > 0. Then V {t ,0) = V {t ,0^) = mo{t) and by dominated convergence for ;<;> 
([T6l ) is transformed into 

V{t,x)=ViO,x)+ ['v'{s,x){-Vis,x) + {moiO)-^is))+Xis))ds (40) 

In the sequel we denote the derivative of functions f{t,x) with respect to the time 
and space variables by f{t,x) and f'{t,x), respectively. 
Let 

U{t,x) :=V(?,x)-(mo(0)-0(0) (41) 
Thus (l40l ) is transformed into 

U{t,x) = U{0,x)+ f -U{s,x)U'is,x) + X{s)U'{s,x)ds + ^{t) (42) 
Jq 

Since V{t,-) is a Laplace transform we have 

U{t,0) = -Q{t) U'{t,0) = -mi{t) limU{x) = -mQ{0)+^{t) (43) 

X — >oo 

and ?7 is a monotone decreasing convex function of the variable x for every t. 

Definition 2.1. Denote by X{t,u) the inverse function ofU {t,x) with respect to x, 
that is U {t,X{t,u)) = u. 

The domain ofX{t,u) in the variable u is {—mo{t) +<t>{t),—Q{t)]. 

X{t,-Q{t))=0 (44) 
The notion of X (?, •) and a version of the following lemma is already present in ||8l. 

Lemma 2.1. IfX{t,u) is defined using a solution of the general frozen percolation 
equation then the following identity holds: 

X{t,u)=X{0,u-<^{t))+t-{u-<^{t))- f\{s)ds+ f'<^{s)ds (45) 

JO Jo 

Proof. We fix an jc,„,„ > 0. For any x > ;c„„„ we have 

|[/(?,x)|<mo(0), \u\t,x)\<'^, \u'\t,x)\<"-^, (46) 

Xmin JC,„/„ 

moreover supo<«7?i(?) < +°°. For an x(0) > x„„„ denote by x{t) the solution of 
the integral equation 

x(t)=x(0)+ [' U(s,x(s))-Hs)ds (47) 
Jo 

This equation is well-posed on the domain x{t) > x,„,„, since U{s,x) — X{s) is 
bounded and Lipschitz-continuous in x. 
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Moreover 

x{t + dt)-x{t) = 0{dt), \U{t,x{t))-U{t,x{t + dt))\ = 0(— ). 

If we differentiate ^ w.r.t. x we get \U'{t + dt,x) - U'{t,x)\ = o{4-). 

U{t + dt,x{t + dt))- U{t,x{t)) - (<!>(? + dt) - 0(0) = 

{U{t + dt,x{t + dt)) -U{t,x{t + dt))) + 
{U{t,x{t + dt)) - U{t,x{t))) - {^{t + dt) - 0(0) = 

rt+dt 

-U{s,x{t + dt))U'{s,x{t + dt)) + X{s)U'{s,x{t + dt))ds+ 

rt+dt Jf2 

U'{t,x{t + dt)) / U{s,x{s))-X{s)ds+0[ — ) = 

/t+dt 
U{s,x{t + dt)) {U'{t,x{t + dt))- U'{s,x{t + dt)))ds+ 

t+dt 



U'{t,x{t + dt)){U{s,x{s)) - U{s,x{t +dt)))ds+ 
r'+'i' , . , . ,df . ,df 



x^. ' 

"mm mm 



/ X{s){U'{s,x{t + dt))-U'{t,x{t + dt)))ds+o{—) = 0( 

Thus U{t,x{t)) = U{0,x{0))+(£>{t), and if we substitute this back into we get 

x{t) =x{0)+tU{OXO))+ f'<^{s)ds- f\{s)ds 

Jo Jo 

By the definition ofX{t,u) we have X{t,U{t,x{t))) =x{t), and by substituting 

u = U{0,x{0)) + ^{t) 
we obtain (|45l). □ 



Since v(0) G V*, V{0,x) is well-defined and analytic for all x£R, thus X(0,m) 
can be analytically extended to (— »io(0),+oo). (|45] ) makes it possible to extend 
X{t,u) to (— mo(0) +0(0,+°°) analytically. The extended X{t,u) is a strictly 
convex function of the u variable. If we differentiate (1451 ) w.r.t. u, we get 



X'{t,u) =X'{0,u-^{t))+t (48) 

Definition 2.2. Define F{t,w) by the identity 

F{t,-X'{t,u)) = -u (49) 

Thus —F{t,w) is the inverse fiinction of —X'{t,u). If X denotes the Legendre- 
transform ofX w.r.t. the variable u, then 

G{t,w) ■.= X{t,-w) = -mm{wu + X{t,u)} =wF{t,w) -X{t,-F{t,w)) (50) 

U 
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Let 

E{t,w) = G"{t,w) =F'{t,w). (51) 

We call E(t, •) tlw critical core of\{t). If we use the extended definition ofX then 
G{t,w) is well-defined and analytic for all w > —t. 

We have 

It follows from the properties of the Legendre-transformation and (l45l) that 

G{t,w) =G{0,w + t)-w^{t)- ['^{s)ds+ [\{s)ds (53) 

Jq Jq 

F{t,w) =F{0,w + t)-^{t) (54) 

E{t,w)=E{0,w + t) (55) 

G{t,-) is strictly convex and G determines X uniquely since the Legendre- 
transformation is invertible. Define 

w*{t) := -X'{t,0) F{t,w*{t))=0 argmin^,G(f,w) =w*(0 (56) 

X{t,0)=0 =^ G{t,w*{t))=0 =^ VwG(f,w)>0 (57) 

e(f)=0 =^ w*{t) = —^>0 (58) 

mi[t) 

oo 

x*{t) = mf{x : y Vk{t)e-'''' < +00} = rmnX{t,u) =X{t,-F{t,0)) = -G{t,0) 

(59) 

3 The frozen percolation equations are well-posed 

Lemma 3.1. The alternating equation ( 1191 ) is well-posed. 

Proof. If we are given the sequence of burning times < Ty < T2 < ... the so- 
lution of (fT9l) can be uniquely constructed by using induction on /: if we already 
have the solution on [0,7^''], then we are given mo{TP), so we can uniquely solve 
the sequence of ordinary differential equations (fT9l ) for vi,V2,... on [7^-*, 7^-^ J by 
repeatedly applying the Picard-Lindelof theorem, since the equation for Vk only 
contains vi , . . . , v/t on its right-hand side. □ 

Lemma 3.2. The solution of the integral equations (116b is unique for every initial 
condition v(0) G V* ifX{t) is nonnegative and continuous. 

Remark 3. Choosing X{t) = implies the uniqueness of the solutions of ([TT] ). 
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Proof. In order to prove the uniqueness of the solution of (fT6] ). we only have to 
prove that given two solutions with the same initial condition, the function <I>(f) = 
mo(0) — mo(f) determined by the two solutions is the same, because mo(f) and ([T6l) 
determines v^.(f) for all k uniquely. For a solution v(f) of ([T6l ) we can define U by 
dHJ, then X by Definition [211. which satisfies (US) and the G of Definition l2!2l 
satisfies (l53l) . 

Assume that G\ and G2 are obtained this way from two solutions of (fT6l) with 
the same initial condition G(0, w). Let G = G\ — G2 and <i> = — O2. Then 



= -wO(f)- f ^{s)ds 
Jo 



Now by Owe have 6(0 = 0, thus (HI) =^ X(f,0) = 0, and ([57]) =^ min„,Gi(f,w) ^ 
minvi,G2(f,w) = and (|58l) =^ w- (?) := aigmm^,Gi{t ,w) > for / = 1,2, thus 
wehaveG{t,w\{t)) < and (5(f,W2(0) > 0. Thus<i>(?) and /g 0(5) cannot have 
the same sign. But if (fi,?2) is a maximal interval such that for ti < t < t2 v/e have 
jQ^{s)ds > then /(J' <i)(5)j5 = and 

?G[?i,?2]^ /" ^{s)ds>0 ^ ^t) <0 ^ [ ^s)ds<0 
Jo Jti 

which contradicts the definition of ?i and ?2- Thus jQ^{s)ds < for all t and 
interchanging the roles of Gi and G2 we get /q <i>(5)(i5 = 0, so (?) = <I>2(?). □ 

Lemma 3.3. If we find a function (p(f) such that defining ^{t) := fQ^^{s)ds and 
G(f,w) ftj (l53l) we have 

minG(f,w) = and w*{t)= argmin,,,G(?, w) > (60) 

w 

for all t, then the solution of (|12l) with the same ?i(-), O(-) and initial condition 
satisfies (fT6l ). 

Proof. Since the Legendre-transformation is invertible, from (l60l ) we get 

X{t,0)=0 and X'(?,0)<0. 

X(f,M) is strictly decreasing for u <0, thus it is the inverse function of an U{t,x) 
satisfying U{t,0) = 0. If we plug <!>(•) into ^ then we get 6(0 = -U{t,0) = 0, 
therefore ([T5] ) is satisfied. □ 



Lemma 3.4. The o/f/je unique solution of (1171 ) w 

ift< 
F{0,T) ift>T 



where 



mi(0)" 



J^^{t)dt = !^ 



//T < r« 

G(0,r) ifT>TS 



(62) 
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Proof. The solution is unique according to Lemma 1X21 and to prove its existence 
we only have to find a function q)(f) that satisfies the criteria of Lemma [33] (with 
X(0 = 0). We will show that 

^{t)=l[t>-^]E{0,t) (63) 

does the job. For t < this is trivial by looking at ^3^: G{t,w*(t)) = and 
^*(0 = 4o)-^>Oif^(0^0. 

We will show that for t>T«we have G{t,0) = and F{t,0) = 0, that is w* (t) = 
0. F(0,r«) = G(0,r«) = by (EH) and w*(0) = = T«. F{t,0) = follows 
from (|54]) and 



^{t)= [\{s)ds= t E{0,s)ds = F{0,t)-F{0,T^)=F{0,t) 
Jo Jts 

By (|53] ) we have 

G(?,0) =G(0,0- ['^{s)ds= [' F{0,s)ds- f F{0,s)ds = 
Jo Jts Jts 



□ 



The well-posedness of the integral equation ([TT] ) implies that of the correspond- 
ing differential equation, since mo{0) — = mo{t) is a continuous function of t, 
thus Vk{t) are differentiable. 

We have shown that the solution of ([17] ) has infinite first moment after the 
gelation time: = m\ [t) = +oo for all t > T^. 

Definition 3.1. Let E{0,w) denote the critical core of\{0) (see Definition \2.2\l . 

<w\<W2 define 

Einf{wi,W2) := min £'(0,w) and £'i«n(wi,W2) := max £'(0,w). 

Wi<W<W2 Wl<W<W2 

Lemma 3.5. Ifw> ^^^^q^ then the inequalities 

"^('^ ' <£(0,w)<i, (64) 



hold. Thus E sup < nii(O)^ and Einf{w) > ^^iM_L 



mi (0) 1 
m2(0) w^' 



For all w > — ^ we have 

— mi(0) 

|£''(0,w)| < 4m2(0)2m3(0) =: D (65) 

which implies 

Esup{wi,W2) - Ei„f{wi,W2) <D-{W2-Wi) (66) 
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Remark 4. If m\{Q) = ni2(0) then the upper and lower bounds in (1641 ) coincide. 
This can only happen if Vk{0) = wii(O) • = 1], this is the case known as the 
monodisperse initial condition (the initial graph has no edges). 

Proof Let U{x) := U{0,x). Recalling ^ E (o,-tj^) = ^-^§M^ holds. The 



U'(x) J ~ U"(x) 

upper bound of dMl) follows from -U'{x) < U"{x), and > U"{x) 

holds because log {—U'{x)) is a convex function, thus > jp^ = ^^^0) • 
bound on the Lipschitz constant (|65] ) follows from 



{UyU'" (Ur /|,..^2..«| , .|,.,.^3 



{u"y u" 



< \{U'yu"'\+3\{Uy\<4m2{0yrm{0) 

□ 



Now we turn our attention to the subcritical equation (fT6l ). We assume > 
for all t. If we substitute x = into the differential equation (l42l ) and assume 
\U'{t,0)\ < +00 then (formally) we get 

4>(0 = cp(0 = -U\t,0).X{t) = m,{t)X{t^ - ^^'^ 



Definition 3.2. If\(0) G V* and X{t) is a positive continuous function then the 
subcritical control differential equation for w*{t) is 

w*{t)= ^ , — - 1 (67) 

w*{t)E{0,t + w*{t)) 

with initial condition w* (0) = — ttvt = T^. 

^ ' mi(0) 

Lemma 3.6. The subcritical control differential equation is well-posed and the 
function 

w* [t) 

{where w*{t) is the solution of (l67l ) with w*(0) = ^^^^j satisfies the criteria of 
Lemma \33\ which implies the existence of solutions to (I16I ). 

Proof. We prove the statement of the lemma on [0, T]. The Picard-Lindelof theo- 
rem and the Lipschitz-continouity property (1651 ) gurantee the existence and unique- 
ness of the solution of (|67] ) before the graph of the solution exits 

{(f,w*):0<?<r, w,;,,<w*<wL„ w*+t>w*m (68) 

for some < w*,,.,, < w*„,,^ < +oo. 

Let Xi„f := info<r<r?i(f), Xvup := supQ<j<7-?i(f). From (l67l) and a "forbidden 

region"-type argument we get that w*(?) +f > w*(0) and w*(?) > min{£^^,w*(0)}, 
since 

w*{t)>0 =^ l(w*(0 + 0>0 
dt 
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w*{t)+t>w*{0) =^ E{0,t + w*{t))<E,,p, 

thus w*{t) < ^ vv*(f) > 0. 

Now we prove that w*{t) cannot grow too fast using the lower bound of (l64l) . 
w*(0 < y{t) where ^(0) = w*(0) = and 



since y{t) is increasing. Thus < a •y{t) +b for some a and depending only 
on the initial data, the function X{t) and T. Thus 

w*{t)<w*{0)e"' + --{e'" -\). 
a 



Now we can see that the graph of the solution of (1671 ) indeed doesn't exit 
until ? = r if we define 

w:,„=min{|^,r^} and wL,:=(r^ + ^K^ (69) 

'-'Slip " 

Now we prove that (p(f) := satisfies the criteria of Lemma [33] by showing 
that 

G(?,w*(0) = and F{t,w* {t)) =0. 

This holds for f = 0, so it suffices to check jjG{t,w*{t))=0 and f^F{t,w*{t)) = 0. 
Using dm) 

at \ w*{t)E{0,t + w*{t)) J w*{t) 

If we combine F{t,w*{t)) = with ([541) we get 

F{0,t + w*{t))=^it) (70) 
It is straightforward to verify jjG{t,w* it)) = by using and dTO]). □ 
This completes the proof of the well-posedness of (fT6l ). 



4 Proof of Theorem O 

We consider the sequence Fpj of probability measures on the compact space Ww[0, T]. 
From Prokhorov's theorem it follows that any subsequence of the measures con- 
tains a sub-subsequence that converges weakly to a limiting measure on Wy^jO, T]. 

Lemma 4.1. Any weak limit point of the measures Pa? is concentrated on the set of 
solutions of the general frozen percolation equation (fT2l) . 
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• If/j{N) = 1, then the X{t) rate function of (|12l) is equal to the X{t) of 

• If^{N) <^ 1, then the X{t) rate function of (|12l) is equal to 0. 
Proof. From (O and dH) it follows that 

Lvf (0 := Jim E(vf + ^0 - vf (0 = 



2 V A^/ l/'^^A^ \ N 

- ^ ((mo(0) - 0^(0) + MOa'(A^)) -vf + vfvL + ^ - 1[2|^] 

(71) 

M(f) = (0 - (0) - jQLv'l{s)ds is a martingale and 

LM^it):= \mi'EiM^it+dt)-M^(t)\!Ft) = lim E((vf (f + J?) k?) < 



N J \\ 2 J N ^ ' ' ' 'J \N 

Thus E (M(r)2) = E {jQLM^{s)ds) = O (^) if we fix k. It follows from Doob's 
maximal inequality that for all £ > 0, A: > 1 and T < +oo we have 



lim P ( sup 

vO<?<r 



vnO-v^O)- Lvns)ds 



> £ = (72) 



If we rewrite this equation in terms of the functions (vv^('))^-i the claim of the 
lemma follows. □ 

Lemma 4.2. Ifjj<^ iu{N), < >i/„/ < l{t) and\{0) G V*, then for any weak limit 
point P of the sequence of probability measures P^r on WyyjO, T] we have 

p(e(0=o) = i (73) 



The subcritical and critical parts of Theorem 11.21 follow from Lemma 14.11 
and Lemma |4^ : any weak limit point P of the sequence is concentrated on 
the set of frozen percolation evolutions satisfying (fT2l) & (ITSl) . When iii{N) = 1, 
P is concentrated on the unique solution of ([T6l ). when ^ <C f^{N) <^ 1 then P is 
concentrated on the solution of (fTTl) . 

In the rest of this section we discuss the proof of Lemma l4!2l 
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Lemma 4.3. We consider a solution of the general frozen percolation equation 
(fnt with initial condition v(0) G V*. If X{t) = or < < X{t) < Xsup < +°° 
then there is a constant C* such that for all t\ < t2 we have 

e(f2)-e(fi)<c*-(f2-?i) (74) 

Proof. First we prove that there exists a constant C depending only on the initial 
data v(0) and ?i,>,y such that 

mi{t)<C (75) 
\fV{t,x) = ZT^y Vk{t)e-'^' then by ^ we get 

V{t,x)= V'{t,x)-{{mo{0)-^{t)) + X{t)-V{t,x)) (76) 

V'{t,x)= V"{t,x){mo{0)-^{t)-X{t)-V{t,x))-V'{t,x)^ (77) 

Substituting V{t,x) - (wio(O) - 0(0) < and ^^^T" - ^"(^■^) ^ we get 

I (-v'(M)) < • (i - ^ (-nM))) 

which implies -V'{t,x) <max{mi(0), =:Cfor allx> and f by a "forbidden 
region"-argument. Thus by letting ^ 0+ we get dTSl ). 
Now we show that for some constant C2 we have 

(V(f,x)-(mo(0)-<I>(0)^^'(^^) <C2 (78) 

for all x>0.1f Xinf <X{t), then by ^ and -mo (0) <V{t,x)- (mo (0) - 0(0 ) < 
we get (|78]l with C2 = mo(0)C. 

Denote by U{t,x) := V{t,x) - (mo(0) - 0(0)- If Ht) = then by ^ and 
(177] ) we get 

^ ([/(?,x)y'(?,x)) = -2V'(?,x)2[/(r,x) - U{t,x)^V"{t,x)+V'{t,x)j^^{t) < 

{-U{t,x))V'{t,x)^ (2--^U{t,x)V'{t,x)) 

Thus we have (1781) with C2 = max{mi(0),2£'i.„p} again by a "forbidden region"- 
argument. Substituting the bounds (1751 ) and (1781 ) into (1761) we get 

j^{-V{t,x))<C2 + C-K,p=:C* 

for all X. Thus y(?i,A;) — V{t2,x) <C* ■ {t2 — ti). Letting x ^0+ and substituting 
into ([T3l ) the claim of the lemma follows. □ 
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We are going to prove Lemma l42l by contradiction: in Lemma \4A\ we show 
that if 9(-) ^ in the limit, then there is a positive time interval such that 9(f) has 
a positive lower bound, and that this implies that even in the convergent sequence 
of finite-volume models, a lot of mass is contained in arbitrarily big components 
on this interval. Than in subsequent Lemmas we prove that these big components 
indeed bum, which produces such a big increase in the value of the burnt mass <!>(•) 
that is in contradiction with <!>(•) < wio(O). 

For any frozen percolation evolution obtained from a frozen percolation Markov 
process on a finite number of vertices we obviously have Q^{t) = (see Q and 
(fT3])). thus 

V^GN £vf(0 = mo(0-<(0='"o(0)-0^(f)-w^(0 (79) 

k>K 

Lemma 4.4. IfFfj^F where P does not satisfy (173] ) on [0, T], then there exist 
El, 82, £3 > and a deterministic t* G such that for every K < +00, every 

m < +00 and every sequence 

f* - £1 < ai < Pi < a2 < P2 < • • • < a„, < (3m < f* 
there exists an Nq < +°° such that for every N > A'o and 1 < i < tn we have 

¥n( max V vf(0>£2 | >£3. (80) 

Proof. First we prove that if P does not satisfy (1731 ) then there exist £i,£2,£3 > 
and £1 < ?* < r such that 

P( inf e(?)>£2)>£3. (81) 

Since (1731) is violated, we have P( supo<,<7- 6(0 > £) > £ for some £ > 0. 

Let L := [^^J and := ^ for 1 < / < L where C* is the constant in (1741) . 

By Lemma 14.11 the random frozen percolation evolution obtained as a weak 
limit point satisfies (fT2l) with a possibly random control function so (1741 ) holds 
P-almost surely for the random element of iV^ [0, T] obtained as a weak limit point. 

Since 0(0) = we have 

{ sup e(0>£}cU{e(?,)>-} 

0<f<r ^ 

almost surely with respect to P. Thus P(6(?*) > |) > f for some t* G . . 
Using (ITU) again ([H]) follows with £1 := £2 := |, £3 = f. 

Now given K and the intervals [a;,P,], 1 < / < m we define the continuous 
functional f : Wjfi, T] ^ R by 

fii{M-)):=, M-)) ■■= r {mo{0)-WK{t)-^{t))dt 

Pi — CLj JUi 
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Thus for all / 

Hi := {((w,(-))r=i M-)) e 'H^w[o,r] : /,-((w,(.))r=i M')) > £2} 

is an open subset of 'FP'wlO, T] with respect to the topology of Definition 1 1.7 1 Thus 
by the definition of weak convergence of probability measures we have 



lim FN(Hi) > F(Hi) > P inf 6(0 > £2 > £3 

W^oo \t*-ei<t<t* J 

from which the claim of the lemma easily follows by dTOl ). □ 

Lemma 4.5. If <^ ^{N) and < ?i,„j < X{t), then for every £2 > there is a 
£4 > such that for every ? > there is a K and an N\ such that 

yN>Ni £ (0) > £2 ^ (<I>^(f)) > £4 (82) 

k>K 

The proof of Lemma l431 will follow as a consequence of the Lemmas [4.61 and 

EH 



Proof of Lemma \4?2\ We are going to show that if there is a sequence Pa? such that 
the weak limit point P violates (173] ) then for some A'^ we have 

(a>^(r)) > mo(0) (83) 

which is in contradiction with (fT3l) . 

We define £1, £2, £3 > and t* using Lemma l4T4l Next, we define £4 using this 
£2 and Lemma 1431 Given these, we choose the. so small that 

^ £3£4 > mo(0). 

We choose K and N\ big enough so that (l82l ) holds for this t. Further on, we fix 
the intervals [a,, |3/], \ <i <m = so that a;+i — P,- > f holds for all / and also 
r - P,„ > r holds. We choose A^o such that ([8OI) holds and let := max{A/o,A^i}. 
Finally, we define the stopping times Xi,X2, . . . ,Xm by 

X,- := P; A min{? :?>«,• and vf {t) > £2}. 

k>K 

We have x,- + ?* < + < a,+i < X;+i. 

Using the strong Markov property, (l82l) and (l80l) . the inequality (1831 ) follows: 



E(<D^(r)) >£E(O^(x,- + f*)-<I>^(x,0) > 

(=1 

m I 

£E ( E((<I.^(x,- + f*) - <I>^(x,))I[i: vf (xO > £2] I :F.,) 

/t>A: / 

££4? ( £ vl (x,) > £2 ) > ni£4£3 > mo(0). 



> 



(=1 \k>K 



□ 
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For a frozen percolation evolution defined by ^ we have 

U{t,x) = £ vf (0^-^- (mo(0) -0^(0) =V{t,x)-m^{t) = £ (0 (.-^^ - 1 



k>l k>l 

(84) 

We will make use of the following generating function estimates in the proof 
of Lemma 1431 

If U{x) = i:i,>i Vk (e"'^-^' - l) where v G V then 



£v^>e =^ Uil/K)<ie-'-l)e (85) 

k>K 

U{l/K) < -£ =^ Y^n> e/2. (86) 



Lemma 4.6. There are constants C\ < +°°, C2 > 0, C3 > such that if 

£vf(0)> 82 (87) 



k>K 



for all N then 



lim P 



( I vf(7) + <I>^(7)>C2S2) =1 (88) 

\k>C3EiN^/^ J 



Where t = ■§^. 

Sketch proof. If we let A'^ — > 00 immediately, then by Lemma 14.11 we get that the 
limiting functions vi (?),V2 (?),.. . solve (fT4l) with initial condition v(0), a possibly 
random control function <!>(?) and some nonnegative rate function 
The ^ 00 limit of (HI) is 

e(z")+0(7) >C2e2 (89) 



Now we prove that if v(-) is a solution of (fT4l ) then Y.k>K'^k{'^) > £2 implies 
with Ci = 4 and C2 = |. This proof will also serve as an outline of the proof of 
Lemma 1431 

In order to prove (l89l) define by (l39l ). Thus V(f,x) solves 

V{t,x) = V'{t,x) ■ (mo(0) - ^{t) + X{t) - V{t,x)) (90) 
Define U{t,x) by (|4TI) . Define the characteristic curve x{-)hy 

x{t)=v{t,x{t))-{m{0)-^{t) + Ht)) 4^) = J (91) 

Let v(0 := V{t,x{t)). Now by ([901) and ^ we get 

v(0 = V{t,xit))+V'it,x{t)Ht) = (92) 
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Thus v(?) = v(0), moreover by dHJ we get U{t,x{t)) - U{0,x{0)) = ^{t) and by 
V{t,x{t)) = V{0,x{0)), y (0,;c(0)) - mo(0) = U{0,x{0)) and ^ we get 

^W = 4+ f^{')ds- f'k{s)ds + t-U{0,]-) (93) 
K Jo Jo K 

By ([85]l we have U{0,j) < -582. In order to prove that Q{t) +(£>{!) > ^82 with 
f = we consider two cases: 

If <t>{t) > ^82 then we ai-e done. If <I>(7) < ^82 define x := min{f : x{t) = 0}. 
By (|93]) we have 

1 - - / 1 \ 1 12^ 



Thus T<t. 



-e(T) = u{x,o) = t/(T,x(x)) = u{o, ^) +o(t) < -^82 + ^82 = -^82 



<I>(f) is increasing. □ 



Thus ^82 < e(T) < e(x) + <t>{x) < e (f) + <J> (0 because by (US the function Q{t) + 



To make this proof work for Lemma 146] we have to deal with the fluctuations 
caused by randomness and combinatorial eiTor terms. 

Proof. Given a frozen percolation evolution obtained from a Markov process by 
© define U and V by dill. 

Using (TtTI) a straightforward calculation shows that 

LV(t,x) := lim ]-E(V(t + h,x) -V(t,x)\ J^t) = 

h^0+ h ^ ^ ' \ / 

V'{t,x) ((mo(0) - 0^(0) + mm - Vit,x)) + 1 {V"{t,x) - V"{t, 2x)) (94) 

Given the random function V{t,x) we define the random characteristic curve x{t) 
similarly to (1911 : 

i(0 = V(?,x(0)- ((mo(0)-<I>^(0) + MO^(A^)) , m = p (95) 

This ODE is well-defined although V{t,x) is not continuous in t, but almost surely 
it is a step function with finitely many steps which is a sufficient condition to have 
well-posedness for the solution of ( |95l ). Define v(f) := V{t,x{t)). 



1+ f {v{s)-v{0))ds+ f'^^is)ds-m f Ks)ds + t-U{0,\^) 
K Jo Jo Jo 

(96) 
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Putting together (l94ll and (1951 ) we get 

^lim lE(v(f + /i) -v(0) I ^0 = ^ (^"(^^(0) - V"{tMt))) = O (^l^V"{t,x{t)) 

(97) 

Thus v(0 = v(f) - /o ^ - V"{s,2x{s))) ds is a martingale and by © 

and ^ we get 

lim |E(v(f + /j)2-v(0^|iFr) = hm \^iiV {t + h,x{t)) -V {t,x{t))f \ft) 



/=i 

Define the stopping time 

Xn := min{f : x(f) = N^'^/^}. 

(Note that we could replace N^^^^ by A'^^''', < y < 1/2 without changing the 
proof.) 

It follows from (l46l ). (|97] ). (1981 ) and Doob's maximal inequality that 



sup |v(f AXa? AT) -v(0)| ^0 as //^oo (99) 

0<t<T 

By ([85]l and ^ we have 

[/(0,x(0))<(e"'-l)£2=:-e5 (100) 

Let 

Ayv:={ |v(^)-v(0)|^/^<-}n{|v(x^Ar)-v(0)|<S5/3}, 

Bn:= {<I>^(x;v)<e5/3}. 
-— 3 3 

^""7i:|[/(o,x(o))| 

We are going to show that that there are constants C2 , C3 < +00 such that 

An<^\ I vf(7)+<I>^(7)>C2S2| (101) 
which together with (l99l) implies hm^r^c^P (A/v) = 1 and 
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First we show that 

AnHBn ^{lN <t}. (102) 

If we assume indirectly that A^v, Bj^ and > ? hold then /q \v{s) —v{0)\ds < ^, 
so by ( |96l ) we get 

x{t) < ^ + f^{s)ds + tU{0,x{0)) < - ^ + ?^ <0. 
K K Jo A3 

But jc(7) < is in contradiction with Xn > t, thus (11021) holds. Assuming Aj^ and 
En we obtain 

|[/(x^,^(x^)) -[/(0,^(0))| < |v(tw) -v(0)| +<I>^(x^) < £5/3 + 85/3 

which together with OOOl) implies AivnBiv ^ {?7(Xa?,A^"^/^) < -£5/3} 
By i 



A/v^(AwnBw)UB^c{ £ vfM>e5/6}u{<I>^M>£5/3} 

/t>A'l/3E5/6 

^1 I vf(x^)+<I>^(XA,)>C2S2[ 

U>C3E2/V'/3 J 

with C2 = C3 = (1 — e^')/6. But Li-x^jEjAfi/s (0 +*J*'^(0 is ^ monotone increas- 
ing function of t, from which (IIOII ) follows. □ 

Lemma 4.7. There are constants C4 < +°°, C5 > such that if 

I vf(0)>C2S2/2 

/or £?// A'^ then with 

In := C4e2 2(Af-'/3log(Af) + iN^liN))-') (103) 

we /lave 

lim E (^^{In)) > C5£2. (104) 

Remark 5. The upper bound (11031 ) is technical: on one hand it is not optimal, on 
the other hand, for the proof of Lemma l?!5l we only need tf^ \ as N ^ 0°. 



Proof. If V is a vertex of the graph G{N,t ) let C/v(v,?) denote the connected com- 
ponent of V at time t. Denote by X;,(v) the freezing/burning time of v. 

^N{t) := {v : |Civ(v,0)| > C3e2A^5 and Xfo(v) > t} 
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We fix a vertex v G ^n{0). 

CN{t) ■= ^\CNiv,t)\ 
^N{t) ■=^\9{N{t)\ 

ZN{t) := ■^|:%(0)\:%(0| =H'w(O)-w/v(0 

Thus CN{t) is an increasing process until Tfo(v), WN{t) is decreasing, ZN{t) is in- 
creasing. We consider the right-continuous version of the processes CA?(f), wa?(?),za?(?). 

Wiv(O) > C2Z2/2 =: £6. 

We are going to prove that there are constants C4 < +00, C5 > such that 

limE(z^(7A,))>C582 (105) 

with In defined as in (11031) . This imphes (11041) . 
We define the stopping times 



: min{f : Wiv(f) < £6/2} 
: min{? : ca?(?) > £6/4} 



Let := C3£2A^3 . Since v G :?/„(0) we have 

Civ(0 > CA?(0) = > - 

If C"/v(v,f) is connected to a vertex in ^N{t) by a new edge at time t then 

/ N \ loef2W 
CNit+)-CNit-)>-, log(cjv(?+))-log(c;v(f^))>log 1+ .J > ^ 



Af' NcN{t-) J ~ NcN{t- 



lim -J-E(log(cw(? + -log(cw(0) I ^0 > 



l0g(2W 1„, , .^ A", N 

• ^ |Civ(v,OI - \cdv,t)\) I{r <x,(v)} > 

log{2)N ■{wN{t)-CN{t))l{t <T,(v)} > 

log(2)A^^I{,<,} = AT V3i^ . C2 . C3 • (£2)' • I{,<x} =: a • I{,<x} 
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Thus \og{cN{t)) — a • (? A x) is a submartingale. Using the optional samphng theo- 
rem we get 

log (mo (0) ) - • E (x) > E (log (c/v (x) ) ) - • E (x) > log (cw (0) ) > - log (A^) 
By Markov's inequality we obtain that for some constant C < +00 

P (x < CA^"^/^£2 ^log(A^)) > ^ 



if N is sufficiently large 

If Zg < Zbiv), then Cw(v,Xg) > ^A^, so E (Xi(v) - x^) < {N^{N)Xi„f)-'^ ^, which 



imphes 



p(T,,AT:h{v)<CN-^/h2^log{N)+C'{Nfi{N))-^£2^^ > ^. 



for some constant C'. Define t of (11031 ) with C4 := max{C,C'}. Using the linearity 
of expectation we get 

E {z(t)) = E ( i £ I{T,(u..)<7} ) > eaP (x/,(v) < t) . 

The inequality I{T„<f}^ < z{t) follows from the definition of x^. 

7 < P (Xv. A Xfc(v) < 7) < P (x,, <t)+P (Zbiv) <t)<E (zit)) - + E (z(F)) - 

From which (11051 ) follows. □ 
Lemma 1431 is a straightforward consequence of Lemma l4.61 and Lemma 14771 



5 Properties of the solutions of the frozen percolation equa- 
tions 

Proof of Theorem [731 It is clear from (l63l ) and (l65l ) that (p(f) is continuous. In 
order to prove (l23l ) we need Example (c) of Theorem 4. of chapter XIII.5 of ||4l. 
By ([55]) 

Y"C/ 0"! = ^ = ^ = ^ 
^ ' ' E{t,0) E{0,t) (p(0 



X(.,«) = ^«^ + 0(.^), Urn-^ = Vm) 

By the Tauberian theorem for any t >T^ each of the relations 

~ 1 

-U{t,x)r^x'-'/^y^ and £v,(0~-T-^i/'-iy2^K0 
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implies the other, that is for any t >T^ 

-U(t,x) I — -- 1 ~ , , /2cp(0 

□ 

In order to compare the solutions of ([T9l ) and ([TT] ) we apply the transformations 

yit)^lJ[t,x)^X(t,u)^G{t,w) (106) 

to the solutions of the alternating equations: 
The integral equation 

lJ(t,x)=lJ{Q,x)^ f -U{s,x)U\s,x)ds + ^{t) (107) 
Jo 

holds, but is constant between burning times and jumps by 6(7;/') at 7;.^ which 
means that the giant component is burnt: 

iim-[/(7;.^+£,o) = iime(7;.^+£) = e(7;.t) =0 

E^O E^O 

By Lemma im the formulae (|45] ). (|53] ). (l54l ) and (|55T ) are vaUd (with rate function 

X(0 = 0). 

In between the burning times tP <t < T^^j we have 

X{t,u)=X{Tl'_^,u) + {t-Tf')u and w) = G(7;.^,w + (? - 7;'')). 

If t-Tl' > w* {tI\) then v(f) is supercritical: 

X'(?,0)>0, 6(0 >0, X(?,-6(0)=0, X'(f,-6(0) <0. 

min^, G(f,w) = still holds, but argminj,,G(f, w) = w*{t) < in the supercritical 
phase. Thus — X'(f,0) =w*{t) is well-defined for all f > for the solutions of the 
equations (fT6l ). (fTTl) and ( fT9l ) as well, moreover ( [59l) holds. For the solutions of 
([T9l ) w*(?) is left-continuous. 

By G(?,w*(0) = 0, (US and dTO]) we get 

f ^{s)ds = G{0,t + w*(t))-w*{t)F{0,t + w*{t)) (108) 
Jo 

for the solutions of ([T9l ). 

If v(0 is the solution of ([TT]) or ^ started from v(0) G V*, then ^ 
holds: the evolution of the critical core does not depend on the rate of lightnings. 
One extra parameter is needed to determine v(f) and Q{t): if we know w* (t), then 

pw rw 

F{t,w)= E{0,t + y)dy and G{t,w)= {w-y)E{0,t+y)dy (109) 

Jw*{t) Jw*{t) 

has all the information about \{t) and Q{t), since the transformations (11061 ) are 
invertible (using analytic extensions). 
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Proof of ClaimUl First assume mo(0) = 1. As a consequence of RemarklH we can 
see that 

1 1 w 

E{t,w)=E{0,t + w) = = ^E{\,j), (110) 

but this is the critical core of and together with w*{t) = for ? > = 

1 the identity Vk{t) = jVk{l) follows. We get the explicit formula for Vk{l) in 
the following way: since X{l,u) = X{0,u) + u, the inverse function of V(l,x) is 
— log(v) + v— 1, thus 

^ ^ k=\ 

where W is the Lambert W function, the inverse function of z i— > ze~. 

If mo(0) / 1 but we still have a monodisperse initial condition then (II 101) still 
holds and for t > -^^^ = we have w*{t) = thus Vk{t) = j^e^'' must hold. 

□ 

Proof of Theorem\Tj\ Let //(w) := F(0,w) -mo(0), thus // (^--^77(1)^^ =-V[{),x) 
by (l52l). Using Lemma [3^ and ([541) we get 

F{t,w)=H{t + w)-H{t) and mQ{t) = F{t,+oo) = -H{t) 

for t > T^. vi(0) > implies limv^eo — = 1, so limj^ocf -//(fw) = from 
which limf^oo?mo(?) = 1 follows. Moreover 

1 — = limt-{H{t-{w+l))-H{t)) = limt-F{t,tw) = limf (f,w) 

W + 1 t^oo f^oo f^oo 

where Vk{t) = tvk{t). This implies the pointwise convergence of the monotone 
functions X'{t,u), X{t,u), U{t,x) and V{t,x) to the desired limit as f ^ oo. The 
convergence of vii{t) to follows from the continuity theorem of Laplace 

transforms. □ 

Proof of Theorem It is easy to check that if v^(f) = Vk{t)e^^*'^*\ then V{t,x) = 
V{t,x + x*{t)), sof*(f) = and w*{t) =0, hut E{t,w) =E{0,t + w) =E{t,w), so 
v(f) is identical to the solution of (fT/l) at time t. □ 

Proof of Theorem \1.10\ If we consider the solution of ([T6] ) with given initial data 
and lightning rate function X{t) >0,0<t<T then (1531 ) provides us with a relation 
between our cost (fQX{t)dt) and reward (Jq <i>{t)dt). 

We prove (1351 ) by considering the cases T >T^ and T <T^ separately. 

According to (l62l ). for T >T^ we get 

o<G-"'^(r,o) = r^"'''{t)dt- r^'"\t)dt+ rx{t)dt 

Jo Jo Jo 

by substituting w = into 
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For T < T«, we want to prove > jQ^'"^{t)dt - J^X{t)dt. Substitute w = 
TS-T into (I53]l. Since G(0, T^) = and (r« - r)<I>™*(r) > we get 

< G'"'\T,T^ -T)<- r ^'"\t)dt + X{t)dt. 

Jq Jo 

The proof of the extremum property (l36l ) is equally simple. □ 

If we want to maximize our cost functional for a fixed T > T^, the optimal 
control is not unique, since the only thing we need for 

r ^'"\t)dt- c x{t)dt= r^"'''{t)dt (111) 

JO JO JO 

to hold is G™^(r,0) = 0: if \{T) is critical at time T, then the value of the func- 
tional is optimal. 

Proof of Remark^ In order to prove (|37] ) first pick an arbitrary X > and solve 
(|67] ) with constant ?i(0 = X. Since w*{t) > and w*(0) = r« there is a < f* < T 
such that w*{t*) = T -t* , and the lightning rate function X{t) =X-I[t <t*] makes 
T a critical time, so (lllll ) holds, thus (|37] ). 

Now we prove (1381 ). By using (11081 ) we have to show that 

G(o, T+w*{T))- (w* (r) - £)f (0, r + w* (r)) > g(o, r) + £F(o, r) 

Using G(0, r + w* (r)) - G(0, r) > w* (r)f (0, r) it is easy to see that < w* (r) < 
sis sufficient for this to hold. If there is a <t* <T such that -X'{t* ,-Q{t*)) = 
T — t* +£, then burning the giant component at time t* we get — X'(f^,0) = T — 
t* +e and -X'{T,0) = w*{T) = e. If not, then burning at time T yields < 
-X'{T,-Q{T))=w*{T+) <£. □ 



6 Proof of the subcritical limit theorem 



In order to prove Theorem 1 1.81 . we need to know more about the solution of (|67] ) 

Lemma 6.1. Ify{t) is the solution of the differential equation y{t 
initial condition yiO) = and t >T^ + clog(^) then y{t) < 2c. 

Proof. The solution of this differential equation is 

'TS-t 



yit) 



1 + W exp 



11. 



(112) 
□ 



where W is the Lambert W function. Thus W{x) <x and our claim follows. 

Lemma 6.2. Ifw*{t) is the solution of (1671) with constant X{t) =X<\ then there 
exist d\ and d2 which depend only on v(0) and T such that 



T^ + diXlog{h<t<T 



w*{t) 



E{0,t) 



<d2X^ 
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Proof. We have a uniform a priori bound w*{t) < w*,^,^ for all X < 1 depending 
only on the initial data and T by ( [69l ). Thus by Lemma [331 we have 

< Ei„f := Ei„f{T + w:„J <E{0,t + w* {t)) 

and substituting this inequality into (|67] ) we get 



w*{t) < 



W*{t)Einf 



1 



(113) 



Using Lemma [6m we get 



t := + 



'^inf 



log( 



EinfT>^. 



X 



<t<T 



w*{t)<2- 



-inf 



Define 



z{t) ■-- 



w*{t)E{0,t + w*{t)) 
X 



1 



Using (|67] ) we get 



z(0 + 



w*{t) 



For t < t < T we have 



1 < z(?) < 2 



sup 



1 ^ 1 ^m/ 

> 



E{0,t + w*{t) 
E'{0,t + w*{t)) 



Einf' w*{t)-2 X ' 



E{0,t + w*{t) 



< 



D 

Einf 



(114) 



(115) 



with the D of (|65] ). Solving the linear ODE (|114l) and using the inequahties (1115b 
we get 



|z(0l<2^exp 



IE 



inf 



2 X 



it-t)]+X 



2D 

E^/ 

mf 



Thus for f > f + -^Xlog (^) we have \z{t)\ = 0{X), which implies 

X 



w*{t)~ 



E{0,t + w*{t)) 



0{V). 



If we combine this with 



X 



E{0,t + w*{t)) E{0,t) 
the claim of the Lemma follows. 



<X^2 



D 



□ 



From this Xm\{t) — E{0,t) = (px{t) — ^crit{t) = 0{X) follows which proves 
(|3T| |. Now we prove (|33] ) using Laplace transforms: 
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Lemma 6.3. Let Ux{t,x) be the solution of d42D with a fixed initial condition 
U {0,x) obtained from v(0) G V* and X{t) = X. Then for any t > we have 

Urn ^ , , \ ^ = , (116) 

x^o [/((f,0) ^/T+^ 

Proof Fix X > and denote the solution of (021) with X{t) =Xhy U {t,x). For all 
f > we have 

X"{t,u)>-^ =^ X{t,u)>-^u' =^ \U{t,x)\ = 0{^). 



We use the shorthand notation E =E{0,t + w*{t)). 



2 

X'{t,u) = -w*{t)-\ hO(M^), X{t,u) = -uw*{t) + — + 0{u^) 

E IE 



U{t,x)=Ew*{t) -J{Ew*{t)y + 2E{x-0{U{t,xy)) 



Ew* (t) - v/ {Ew* {t)f + 2Ex + O (x) 



-1 / N -1 

+ 0(1)= ; = + 0(1) 



E 

Because of Lemma [6!2l we have 



w*{ty + ^x+0{x) Jw*{tY + ^x 



lim 



x"hE{0,t + wl{t))E{id,t)wUtY 



from which the claim of this lemma follows. 



□ 



The r.h.s. of (11161) is the Laplace transform of the r(i, 1) distribution and the 
r.h.s. of (1331 ) is the distribution function of the r(^, 1) distribution, so (1331 ) follows 
from the continuity theorem of Laplace transforms. 



Proof of Theorem [73] First observe that instead of proving uniform convergence 
of <i>n to ^crit wc Only need to show convergence on [0, T] for any T, because 

T>T^ ^ mo(T) = [ E(0,w)dw < [ -^dw = - 
Jt+w'{T) Jt w^ T 
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by dMll, thus < <I>„(f) -<I>crif (0 < 7^ for ? > r. If we prove that w*{t) is small for 
t > then we are done by dTOl ) and Lemma [331 since 

0<^„{t)-^cri,{t)=F{0,t + w:it))-F{0,t)<w:{t)E,up t>TS (117) 

^n{t)<^n{TS)=F{0,T^ + wl{TS))<wl{TS)E,up t < T« 

We can give an upper bound on w* (t) for ? > if we replace with Xsup in 
(I67]l: using ([TT3]) we get w*{t) = O (?ilog(i)) if we substitute t>T^ and c = ^ 
into (|112l) . thus lim„^ooW*(f) = uniformly for <t <T. 

We obtain lim„^oo v^(f) = Vk{t) for k = 1,2, .. . by the uniform convergence of 
mo(f) and X"{t) to the critical mo{t) = mo(0) - ^{t) and X{t) =0m □ 

7 Proof of the alternating limit theorem 

We turn our attention to the proof of Theorem 1 1.41 and Theorem |1.9l 

In this section we assume mo(0) = 1 but the results generalize easily to the 
mo(0) / 1 case, since if \{t) is the solution of ( fT9l ) with burning times ,T2^. . . 

then mo(0)v(mo(0)f) is also a solution of ( fl9l ) with burning times ^^^^^ , ^^^^.^ , . . . 

1 



nil 



Definition 7.1. -//'v(?) is a solution of fl79l). fe? w;^(f) := 

//w*(0 > then w%{t) = w*{t), but ifw*{t) < then w^it) = -X'{t,-Q{t)). 

If ns a burning time then w*(f+) := lims^^QW* {t + e) =w*^_{t). 

Lemma 7.1. We consider the solution ofl\19i on [0, T] with an arbitrary sequence 
of burning times. IfT^<t<T and w*{t) < then 



//w*(0 <Oandif 



wi + < r < W2 - /l^lEi!^ (120) 

y £,•„/( wi,w2) 



'^4^w*(0 < < - J^H^w*(o (121) 

-2Einf{wuW2)w*{t) < 6(0 < -2£,„p(wi,W2)w*(0 (122) 
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Proof. By (US, w;(?) = -X'{t, -6(0), (fT09l) and ([61 we get 

ReaiTanging (l50l ) and using (l49l ) we get that w = w;^(?) is the positive root of 
the function 

f{w) ■= G{t,w) -F{t,w)w = Git,0) + (-£yE{t,y)] dy = m + (/(w) -/(O)) 



We prove (fTT9l ) by considering the cases < and > 

separately. 



jf < 1^^, then we prove that w;(0 < 2^'^^ \w*{t)\ by showing 
that/(0) < |/(w)-/(0)| withw = 2yg||w*(0|. 

/•O /•|v>'*(OI y 

m = / (-3;)£(0,?+3;)^/3; < / jr^dy 

J Wit) Jo [t-y) 

by (fT09l) and (il- 

'■'^ ^ •'^ ^' -yo m2(0) (f+3;)2 io m2(0) (^ffl +3;)! 
It is straightforward to check that 



, */ M „ \w*(t)\ 1 /mi(0) 3; mi(0) 



which is sufficient for /(O) < \f{2,/^^\w*{t)\) - /(O) 



4\/m2(0) (f-3;)2 - m2(0) (fM+3;)2 

to hold. 



If^>lViI>then 



/(O) = G{t,0) = f F{t,y)dy < \w*{t)\ < T 

Jw*{t) 



since by (l52l ) we have F{t,y) < mo{t) < mo(0) = 1. Calculating the middle integral 
in ([T23]) we get that in order to have /(O) < |/(w) - /(0)| 

^iog(i + ^^)-i)>r 



m2(0) V f 



is sufficient. Rean^anging this and using > | y ^^^^ we obtain (11191) . 

The proof of the upper bound of (|121l) is similar: using (11091 ) we get that wi < 

f — Ivv* (?) I < ? + w < W2 imphes 

/(O) < ^£™;,(W1,W2)W*(0', /(W) -/(O) < -^£,-„/(wi,W2)w2 
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Using (fT20l) the inequality / y- \J^^^^w* {tyj < follows. The lower bound 
of (|121l) is verified similaiiy. 
If [-6(0,0], then 

X{t,u)<-W*{t)u+]-—-^ 

2Einf\Wi,W2) 

since X"{t,u) with u G [— 6(?),0] is equal to ^(^Q\_^_y•J for some 



thus f +3^ G [wi , W2] by (I120I ). This implies the lower bound of (|122l) . and the proof 
of the upper bound is similar. 

□ 

The proof of Theorem 11.41 is similar- that of Theorem 11.3 1 : if £ = sup,{7^-^j — 
tP} and Tj' <t< T.'[^ then w*{t) = w\{tI') -{t- Tj") > -e and by (fTTOl) we have 
w;(7;.^) = 0(|w*(7;.^|)) = 0{e) on [o,r]. 

Lemma 7.2. We consider the solution of ([19\l with initial critical core E{0,w). 
If Tf < T2 are two consecutive gelation times, then the unique burning time in 
between and T2 is 

f^lyEiO,y)dy 

r"(rf,r/) = 4^ (i24) 

ff, E{0,y)dy 



Q{T\TfJi))=rE{0,y)dy (125) 



Moreover 



Proof. needs to satisfy T2 —T^ = w*^{T'^), but by the proof of Lemma ITTTl 
w\{T^) is the unique positive root of G(r^,0) - /o"'y£(0,r^ +y)dy. G{T\0) = 

- f°g_p,yE{0, T'' +y)dy by so 3'^(0, +y)dy = must hold, from 

which (11241) easily follows. 

By dllll, w\{t) = -X'{t, -6(0) and (fT09l ) we get 



Q^j.i,^^j,^j.t^^*^j.t^^^ r^^^E{0,T' + y)dy= E{0,y)dy 

Jw*{T'') JT^ 



□ 



Definition 7.2. If\{t) is the solution of the random alternating equations (see 
Definition \1.6\ ). denote by < < ... the sequence of random burning times 
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and by = Tf < T2 < .. . the sequence of random gelation times. Indeed Tf < 
TI' <Tf <T^<... 

Let Xi := 7^-^ J — Tf be the length of the i-th critical interval. 

N{t) ■= max{/ : Tf < t}, x(?,/) := Xa,(,)+; 

T(f,0) is the length of the critical interval containing t. 



Let Q{t,i) := 9(7]^/,n_|_j), thus Q{t, 1) is the frozen mass of the first giant compo- 



nent bom after t. 



N(t)+i - ^N{t)+i - y^N{t)+iJ 
g Tb 



w\{t,i) :- 7]v(()+,-+i - 7]v(()+,- - w+(r^(f)+,.) 

Definition 7.3. A nonnegative random variable X has Rayleigh distribution with 
parameter a, briefly X R{o), if 

p (X > x) = exp{-:^x^) -.R{a,x) 

E (X) = Oy^. Y has a size-biased Rayleigh distribution with parameter o, briefly 
Yr^R,b{o)if 

, E(X-I\X>y]) , , 
P(^>^)= E(X) =^-^(^'^) 

The scaling identities 

R(o,x) = R{ao,ax) and Rsb{<^,x) = Rsb{a(5,ax) (126) 
are valid for a > 0. 

The r.h.s of ^ is /?,/,( 

The Rayleigh distribution emerges in our setting in the following way: if we 
consider the solution of the random alternating equations with burning times de- 
fined by a homogenous Poisson process with rate X, forget about the error terms in 
(fml) by assuming wi=W2 then %{t)=lE-{t- Tf) if T;.^"^ <t<Tf, so 

Tf - Tf > w) = exp(-X lEsds) = R{-j==,w). 

From 6(7;*) = IE ■ {Tf - Tf) and (fT26]) we get 6(7;*) ~ R{\f^)- Assuming wi = 

W2 in (ImT) we get w\{Tf) = -w*{Tf), thus X; ~ R{\f^)- 

Lemma 7.3. Ifyjt) is the solution of the random altemating equations with con- 
stant rate function X{t) = X then for every T^ <t <T we have 

E (6(r* )I[r* < T]) = 0{X-^-) (127) 



^[%)+i^T-tj = 0{X-2) (128) 
asX^ 0° where the constant in the O depends only on the initial data and T. 



1 , 
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Proof. Let y{t) : = f - T^^^^y Then 

Ji^O ^^^^ + "^^^ " ^^^^ I = ^ - '^(^) Ji^o - ^ ^ + i^r 

i-Y(0e(r^r4),0)X^r^r4),.) 



2 F 

'-'sup 



by Lemma 17.21 Taking the expectation of both sides of the above inequality and 
applying Jensen's inequality we get 



|e(y(0)<i-^x^|^e(y(0)^ 

This differential inequality together with y(7"^) = implies 



by a "forbidden region"-type argument. Now we prove 

E(7;vV,Ar-r4)) = 0(X-i) (129) 

from which (fT28l) trivially follows. We obtain ([127] ) using ( [T29l ) and 6(r^(y)) < 
lEsup ■ {T^(^^^ - by the upper bound of (1122b . 



A 7^ - 4V) A t) i[t < r^)] + (r^) AT-t)i[t< r^)] 



C{T,y{0))w*_it,0)I[t > r^)] < C(r,v(0))Y(0 
where C(r, v(0)) is defined in (fTT9l) . 

c(r,v(0))Y(0 +c(r,v(0)) (r^, a r - A i[? < r^,] 
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By ^ and (11181) we have 

E ( (^4) AT-t)i[t< r^)]) = E ( (r^) A r - ^) V 



E (7^4)+! AT- j = O (E (Y(0 ) ) + o (E ( (r^) A r - f j V 0) ) = o (X- ^ ) 

□ 

Sketch proof of of Theorem \1.9\ Our aim is to make the following argument rigor- 
ous: Let 

If 1 <C ?i then e(f,l),e(f,2),... ,Q{t,n{X)) are "almost" i.i.d. with distribution 



e(.,/)~/?(V^).x(.,/)«fi) 



so 



i=i 



by the weak law of large numbers. Substituting x = 2\J ^^^^ into 

<l>(f + 8(X),i)-<I>(f,x) i:"iy9(f,/)-I[9(f,/) >x] E(9(f,l)I[9(f,l) >x]) 

lSW^O ^ E(9(?,l)) 

we get dMll- □ 
Proof of Theorem 17791 We use the notations of Definitions 17.21 and l7.3l 

£:=£(f,0)=£(0,0=cp„.Y(0 

We fix a; > and define 

x:=2y-x, 9(f,/,f):=9(f,/M9(f,/) >x], := [6(X)y — (1 +z)J 

By the assumption ^ 5(?i) we have \imx-^oon(k,z) = +°° for any — 1 < z. 
Letm(?i) :=A^(f + 6(?i))-A^(0-l. 

+ 5(X),x) -<I>(f,x) = 9(f,0,f)I[r4) > f] + 

m(;L) 

£ 9(?,/,f)+9(? + 5(X),0,x)I[r4^g(5^)) <? + 6(X)] (130) 

;=1 
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In order to prove (l34l ) we only need to show that we have lim;^^„P(B(X,8)) = 
1 for every £ > where 

because the first and the last term on the r.h.s. of (11301 ) divided by Ed{'k) converge 
to in probability asX^ oo by ([l27l)and?i-2 «;5(X). 

E,,,p{X) := E,,p{t,t + 25{X)), EinfiX) := + 25(X)) 

By (l66l ) we have 

E,up{^) <E + 2D^X) and E - 2D^X) < Ei„f{X) (131) 



]imi_^^C"{X) = lim;,^„C'(?i) = 2, since 5(>i) < 1. 

We are going to couple the random variables T^^^^^.^,w*_{t, l),w*_{t,2),... to 

wL(l),w'_(2),... and w'i(l), (2), . . . 

where w'_{i) ~ R{ , J ) are i.i.d. and w"_{i) ~ /?( ,,,, ) are i.i.d., more- 
over the auxiliary random variables are independent from T^,f\^i- If we define the 



events 



A"(^,z,Z2) := jr^)^, w'la) <? + 6(X) -(1+22)1 

'4'(^,^,^2) := 1^4)^1 +C'(X) •'£w'_(7) >t + d{X) ■ (1 +Z2)| 

then it is an easy consequence of (1128b . X^i ^ 5(X), and the weak law of large 
numbers that 

-l<z<Z2 =^ limP(A"(?i,z,Z2)) = 1 

z>Z2>-l =^ limP(A'(?i,z,Z2)) =1 
Our coupling is going to satisfy 

A"{X,z,l) C n {w'_{i)<w*_{t,i)<w"_{i)} (132) 

for any z- 
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The joint construction of wL{j), w*_{t,j) and w"_(j) for 7 = 1,2,... is as fol- 
lows: given T^(^t)+i andH'l(f, I), . . . ,w*_{t,j — I) we can determine T^^^^^^. by solv- 
ing For s > Let 

lj{s) := Xe(r^(^)^^. + iJi{s) ■= X2Es„p(k)s, /Juis) := X2Einf(k)s. 

Let w*_ (tj), wL (j) and w"_ (7) be the horizontal coordinate of the leftmost point be- 
low the graphs of /u, f^i and jUu of the same standard uniform 2-dimensional Poisson 
process on the first quadrant of the plane. Thus w'_{j) ~ /?( , ^ ), w"_{j) ~ 

\/2£n,;,(/.)^ 

j?( , ' ) are independent from everything that was constructed earlier and 

P (wL(;) < w'LU)) = 1- The joint distribution of T^^^^^^,w*_{t, 1), . . . ,w*_{tj) 
agrees with that of the solution of the random alternating equation. 

We are going to prove (11321 ) by induction. Assume that A"(?i,z, 1) holds. If 

f]{w'_{i) < w*_{t,i) < w^/)} n f]Mt,i) < C'{X)-w^Lii)} (133) 

i=l i=l 

holds for some j < n(k,z), then 

1= 1 1= 1 

which implies iu„{s) < ili{s) < iJi{s) for < 5 < w"_{j) by (fT22]) and 1). 
From this wL{j) < w*_{t,j) < h'"_(7') follows, and (11211 ) can be applied to deduce 



Thus we can replace j with 7 + 1 in (11331 ). This completes the proof of (11321) . Let 

e"(f,/,f) := 2£,„p(X)w"(/) •I[2£,„p(X)w"(/) > x] 

e'(f,/,f ) := 2£,-„/ • l\lEinf{'k)W{i) > x] 
(fml ) and (fT32l ) imply 

A"(>^,z, 1) ^ n < e(f,/,f) < e"(?,/,f)} 

(=1 

B (X,z,s).-| ^g^^^^ <RA^,x) + ej 
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The law of large numbers, (11261 ) and (11311 ) imply that 

z<E ^ limP(B"(?i,z,£)) = 1 and - e < z ^ lim P (b'(>i,z,£)) = 1. 
We can use (11321) and (11211) to show 

A"(X,z,l) C f] {c'{X)wL{i) <x{t,i) <C\X)w'Ui),} 
1=1 

Since 

=max{7 : r^(,)+i + It(?,/) < ? + 5(X)} 

and A"(?i,z,0) C A"(?i,z, 1) by definition, 

A"(X,z,0) C {miX) > n{X,z)}, A'{X,z,0) nA"{X,z, 1) C {m(?i) < «(X,z)}, 

A'(X,|,0)nA"(X,|,i)nB"(X,|,s)nB'(X,-|,s)nA"(X,-|,0)cB(X,s) 
This completes the proof oflimx^„P{B(X,e)) = 1. 

□ 
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